An optically resolvable Schrodinger's cat from Rydberg dressed cold atom clouds 
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In Rydberg dressed ultra-cold gases, ground state atoms inherit properties of a weakly admixed 
Rydberg state, such as sensitivity to long-range interactions. We show that through hyperfine-state 
dependent interactions, a pair of atom clouds can evolve into a spin and subsequently into a spatial 
Schrodinger's cat state: The pair, containing 40 atoms in total, is in a coherent superposition of 
two configurations, with cloud locations separated by micrometers. The mesoscopic nature of the 
superposition state can be proven with absorption imaging, while the coherence can be revealed 
though recombination and interference of the split wave packets. 
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Introduction: When and why mesoscopic objects begin to 
behave according to our classical intuition, as exemplified 
by Schrodinger's famous thought-experiment [1 j, remains 
one of the fundamental questions in physics. Experimen- 
tal progress to demonstrate quantum coherence in meso- 
scopic systems is impressive, with recent creation of su- 
perposition states of macroscopic Josephson currents [2], 
ten-qubit photonic cat states [3] , six-qubit atomic hyper- 
fine cats [4] , interference of fullerenes and even large bio- 
molecules [5] , superpositions of photon coherent states [6] 
and many more [3 [8] . 

In most of these experiments the quantum mechani- 
cal superposition does not pertain to an intuitive classi- 
cal observable taking common- sense values, such as the 
original "alive" or "dead" of Schrodinger's cat. Instead, 
the superposition typically is achieved with intrinsically 
quantum mechanical degrees of freedom (hyperfine- or 
photon number states). Realizations of position-space 
superpositions have been limited to small derealiza- 
tion lengths (several lOOnm for Ref. [5]) the resolution 
of which requires sophisticated near-field interferometry. 
Here, we propose a Schrodinger cat superposition in the 
relative distance of two ultra-cold atom clouds more than 
10/im apart. The relative distances of the two superposed 
configurations also differ on a micrometer scale, hence the 
existence of the two possible cloud configurations can be 
revealed with direct absorption imaging and the coher- 
ence of that superposition can be proven by interference 
upon recombination, taking Schrodinger's cat into the 
optically resolvable micrometer domain. 

In contrast to prior proposals with ultracold or Bose- 
Einstein condensed atoms (e.g. [9HT7]) we use Rydberg 
states, taking advantage of their inherently strong long- 
range interactions and short dynamical time-scales [THJ 
[T9] . The resulting internal forces [20j [21] let the system 
turn itself from a spin cat state into a spatial cat state. In 
addition, Rydberg systems typically allow for an accurate 
control of decoherence mechanisms. 

The scheme (see Fig. [T]) is based on a pair of atom 
clouds, each containing about 20 alkali atoms, which can 
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FIG. 1: (color online) Schematic of Rydberg dressed atom 
clouds, (a) All atoms are in either of two hyperfine ground 
states | g ) , \h). Dressing lasers can couple one atom per 
cloud to either of the | s), \p) Rydberg states. The Rydberg 
states participate in state changing dipole-dipole interactions 
\sp) <H> \ps). (b) Two atom clouds with width a, separated 
by a distance r^i indicates the blockade-radius. Due to 
hyperfine state dependent inter-cloud forces, a suitable initial 
state evolves dynamically into a non-classical position space 
superposition state, with the pair of clouds in either the red- 
or green shaded configuration. 



be in one of two hyperfine levels | g ) and | h ) of the 
atomic ground state. To induce long-range interactions 
between the clouds we weakly dress the states | g ) and 
| h ) with Rydberg states | s ) and | p ) , respectively [22]- 
[25] . These are chosen such that each cloud is in the full 
dipole-blockade regime [26] , where only a single Rydberg 
excitation per cloud is possible. However, the inter-cloud 
distance is so large that excitations in different clouds 
do not block each other. Such interactions can lead to 
collective relative motion of the clouds, with a repulsive 
or attractive character depending on the total hyperfine 
state. 

To realize this scheme, we first identify a suitable ef- 
fective state space and Hamiltonian for our system. We 
then show how to create a hyperfine state, formally al- 
ready a spin cat state [27j [28] , in which the two clouds 
evolve as a coherent superposition of attractive and re- 
pulsive dynamics. After a brief dwell time, single-shot 
absorption images would at this stage show either the 
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green or the red configuration in Fig. [T] To see the co- 
herent character of this many-body state via interference 
fringes, recombination of the two configurations is finally 
possible with the help of an external (double- well) poten- 
tial, as demonstrated in the last section. 
Ultra cold Rydberg dressing, dipole-dipole interactions 
and blockade: Consider an assembly of 2N neutral atoms 
of mass M located at positions r n , restricted to one di- 
mension and confined to a double- well atom trap. Half 
the atoms are localized in one of the wells, forming cloud 
A and the rest in the other well, forming cloud B. Near 
the centres of each well at x = ±d/2, the potential is 
approximately harmonic: V(r n ) = Muj 2 (r n ± d/2) 2 /2, 
and the atoms are initially in the Gaussian trap ground 
state of width a = ^H/Muj. We consider four essen- 
tial states in 8T Rb atoms. Two of them are long lived 
hyperfine states | F, mp ), namely \g) = | 1, — 1 ) and 
| ft ) = | 2, 1 ) (F is the total angular momentum and mp 
the associated magnetic quantum number). The other 
two essential states \v,l) are Rydberg states, designated 
by \s) = | 80,0) and \p) = | 80, 1 ) (y is the principal 
quantum number and I the orbital angular momentum). 
The Rydberg states are coupled to the ground states with 
Rabi frequency ft and detuning A, as sketched in Fig. [T] 
The coupling is off-resonant, hence a = Q/2A <C 1. 
As shown in [24 this arrangement gives rise to effective 
long-range (state changing) dipole-dipole interactions of 
the form D(r)(\gh)(hg\+c.c), between dressed ground 
states | g )~| g)+a\ s), | ft )H h)+a\p). We have D(r) « 
a 4 /i 2 /r 3 , where the transition dipole \i parametrizes the 
strength of the bare dipole-dipole interaction. Hence, 
we can further reduce the essential electronic state space 
of a single atom to | g ) and | ft ) , on which we build the 
many-body basis | k) = | k\ . . . k2N ) = | &i )<£)•• *®| &2Ar ), 
where kj G {S} = ft} describes the electronic state of 
the atom j. We formulate the many-body Hamiltonian 

H = H + Hint, fli„t = Yl D nl {K)a { ^af_+h.c, 
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\k n )(k' n \ where k n: k' n G S in i^ int , while 
D n i(H) = D(\r n — ri\) describes the induced transition 
dipole-dipole interactions. Here the operator acts 
only on the Hilbertspace of atom n and as unity oth- 
erwise. The vector R = {n, . . . r2Ar} T contains all atom 
co-ordinates. Note that there are no interactions between 
two atoms in the same cloud, since the required doubly 
Rydberg excited intermediate state is strongly energeti- 
cally suppressed through the dipole blockade [29 . The 
total number Nh of atoms in state | ft ) is used to classify 
the electronic states, since Nh is conserved by ^i n t- 

Having set up our effective state space and 
Hamiltonian, we can construct adiabatic Born- 



Oppenheimer (BO) potential surfaces J7/(R) defined by 
fr int (R)|<#(R)) = Ui(R)\cpi(R)) [3(K As discussed in 
previous work [3TJ [32] the motion of atoms is deter- 
mined by these BO potentials, as long as non-adiabatic 
effects are small, see also supplementary information 
[33] . We characterize Born-Oppenheimer surfaces in 
the vicinity of the initial configuration sketched in 
Fig. [T] with the dichotomic central many-body position 
Ro = (~d/2, ••• ,-d/2,d/2, ••• ,d/2) T around which 
the positions of the atoms (R) are randomly distributed 
with width a. We choose d = 11/im and a = 0.5/im. 

In Fig. [2] (a) we show cuts through BO surfaces for 
states with Nh = N (half the atoms in | ft)) as a func- 
tion of d. The insets show coefficients Ck of the two eigen- 
states I ^ r ep/att ) = c k,re P /att I k ) with the largest ab- 
solute eigenvalues /7 rep / att (R) ~ ±N 2 D(d)/2. These 
states are of particular interest, since on the correspond- 
ing BO surfaces the entire clouds attract or repel, as de- 
duced from the gradient of U rep / att . Consequently, af- 
ter preparing the twin atom clouds in a hyperfine state 
I ^cat ) = (| ^att ) + I ^re P ))/ V% one obtains a spatial su- 
perposition state as sketched in Fig. [l] through motional 
dynamics. 

If our underlying basis is mapped onto a spin system 
[33] . this process can be viewed as conversion of a col- 
lective spin Schrodinger's cat state into a spatial one. 
The states | ^rep/att ) are close to coherent spin states in 
this picture, as sketched in Fig. [2^i. This conversion does 
not require external fields, but proceeds entirely through 
internal interactions within the system. Note that the 
collective cloud motion in a blockade regime crucially re- 
lies on the dressed character of the interaction. For bare 
dipole-dipole interactions only a single atom per cloud 
would be accelerated [29] . 

Having established the fundamental mechanism be- 
hind our cat, we will outline how the initial hyperfine 
state I ^ cat ) can be prepared, and then proceed to model 
spatial dynamics and interference. 

Initial state creation: The first stage of assem- 
bling l^cat), starting from the simple state |g) = 
(| k) with k n = g. Vn), is to create | ^ rep ). This can be 
achieved on time scales shorter than that of atomic mo- 
tion by using a micro-wave field which couples the two 
hyperfine-ground states so that the atom-field interaction 
Hamiltonian during initial state creation is [34] : 



#ini(R) — H[ n t{H) + i7 r f, 



E 



[^)agV2 + h.c.-A r Kt)a^]. (2) 
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When we analyze the spectrum of Eq. (11) for constant 



i^int and Rabi-frequency ^ r f, as a function of micro wave 
detuning A r f , we see that the eigenstate | g ) at large 
negative detuning evolves continuously into | ^ rep ) at 
A r f = 0. This state is adiabatically followed in Fig.[2](c), 
using the chirped micro- wave pulse shown in Fig. ^T(b). 
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FIG. 2: (color online) (a) Born-Oppenheimer surfaces Uk(R) 
as a function of d for the case 2N = 8, a = 1/xm, a = 0.15 
and fi — 5760 atomic units. Positions R have been chosen 
around Ro in accordance with a single realization of a Gaus- 
sian distribution (width a). The insets show the eigenstates 
belonging to the two colored energies, red, top inset: Coef- 
ficients Ck rep = ( k | ^rep ) of the most repulsive state, blue, 
bottom inset: Ck, a tt — (k| ^att ), most attractive. A plot of 
the modulus |ck a tt| coincides with c\^ rep , since the states dif- 
fer only by signs of coefficients (blue dashed, top inset). The 
transparent spheres visualize the Q-functions of | ^ rep / a tt ) m 
a pseudo-spin picture, in which the states resemble coherent 
spin states [33 . (b-d) Creation of | ^rep } for d = 11/im, using 
a chirped micro- wave pulse as described in the text, (b) Time 
dependence of micro wave Rabi frequency fi r f (black) and de- 
tuning A r f/10 (red), (c) Resulting energy spectrum of H[ n [, 
Eq. ( |11| ) the red line is the state to be adiabatically followed, 
(d) (black) Energy gap between the two highest states of (c), 
compared with analytical prediction Eq. (|3| (red-dashed). 



The pulse avoids non-adiabatic transitions since the pulse 
length T r f is long compared to the inverse energy gap 
AE~ X between the two highest energy eigenstates. The 
latter is well approximated by 



AE(t) = v/A rf (t) 2 + tt rf (£) 2 + g{N)D{d), (3) 

with an only weakly TV-dependent factor g(N) ~ 0.5, 
as shown in the supplementary material [33 j. The re- 
sult Eq. (|3| simplifies the determination of realistic pa- 
rameter regimes. We have numerically modelled the 
pulse of Fig. [2] (b) for 2N = 8 and found a fidelity 
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0.85 when averaging over the 



atomic position distribution. Our creation scheme for 
| \I/ re p ) closely follows the method of [35] . 

The second stage of initial state creation is to convert 
| ^ rep ) into | ^ cat ). We find that | ^ rep ) and | tf att ) are 
always related as shown in Fig. [2] (a): | ^ at t ) is obtained 
from | \I/ re p ) by a 7r phase shift to every coefficient of 
basis states involving an odd number NhA, of atoms in 
| ft ) in cloud A. By applying this phase shift conditional 
on some control atom in a (| ) + | 1 ))/v // 2 superposition 
we achieve our goal. This can be realized precisely as 
in a recent proposal for mesoscopic Rydberg quantum 
computation gates [36], see also [33]. When modelling 



this final step of the initial state creation sequence, we 
find that fidelity loss is negligible compared to the one 
incurred in the previous stage of creating | ^ rep ). This 
situation should persist for larger TV [36] . 
Spatial cat state creation and interference: To turn the 
electronic state | ^ cat ) prepared so far into a spatial cat, 
we keep the dressed interactions switched on for an ac- 
celeration period r acc ~ 6/is, after which they are adia- 
batically switched off to avoid spontaneous decay of Ry- 
dberg population. After mechanical evolution in the trap 
for a time r trap /4 = 7r/(2o;), the clouds reach their max- 
imal displacement, where the macroscopic spatial super- 
position character of the quantum state can be shown 
with /im resolution atom detection. An absorption image 
would always show two inert clouds, with 50% probabil- 
ity at either of the two configurations marked / and O in 
Fig. [3] If instead the spatial dynamics is allowed to pro- 
ceed until time r trap /2 where the spatial wave function 
recombines and all atoms are reunited in the same hyper- 
fine state [42 , we form an interference pattern, demon- 
strating the coherence of the superposition. 

We solve the Schrodinger equation as in [3TJ [33l [37] to 
model the quantum dynamics of acceleration, splitting 
and recombination for 27V = 4 in a plane wave basis and 
for 2 N = 6 in a Hermite- Gauss basis. Interference fringes 
develop in the probability distribution p(z) of the rela- 
tive inter-cloud distance z = (^2 neA R n — ^2 n eB Rn)/N. 
We extract p(z) from the many-body wave function as 
p(z,t) = J dR\ty(H, t)| 2 , where J dR denotes integration 
over all coordinates orthogonal to z. At t = T trap /2, we 
find full contrast interference fringes in both cases. We 
thus believe that they persist also for larger atom num- 
bers, as no new physics enters beyond 3 atoms per cloud. 
Atomic densities and interference for 27V = 4 are shown 
in Fig. [3] which additionally includes results obtained 
with Tully's quantum classical algorithm [20j [32j [38], 
with which slightly larger atom numbers can be treated 
(~ 8 — 12). We find that nonadiabatic effects during the 
acceleration phase are negligible, with a population loss 
of P na = 1 — £> rep = 10 -9 out of the target state | ^ cat ) 
for the situation of Fig. [3] For larger N the situation 
improves further. 

While computational demands limit the simulations 
shown to 2 N < 6, 8 we extrapolate that spatial cat states 
are realistic for up to N = 20 for the parameters used in 
this article. Nonadiabatic effects during acceleration and 
initial state creation are under control for larger N. The 
main limitation comes from the life-time of the Rydberg 
states used for the dressing, since just a single decay has 
the potential to destroy the fragile cat state. However, we 
can choose parameters for which the probability of even 
a single decay is small. This is for example achieved for 
d = 11/im, a — 0.5/im, a = 0.15, assuming 87 Rb atoms 
with v = 80. We use \i = \±§v 2 ! , where po = 0.97a.u. for 
87 Rb. The overall life-time of the system under dress- 
ing interactions is rnf e = r/(2iV), with r = To/a 2 and 
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FIG. 3: (color online) Evolution of the initial hyperfine state 
| ^cat ) into a spatial superposition of cloud locations for 
2N — 4 and trap-frequency u — (2tt) x 400Hz. We com- 
pare the total atomic density n(x,t) from quantum-classical 
simulations (grey shading, black 3D lines), to full quantum 
solutions (yellow dashed). The lines overlayed on the grey 
shading are exemplary repulsive (green) and attractive (red) 
quantum-classical trajectories. Transparent spheres show the 
conversion of a spin cat state into a spatial cat state that 
would under ly this process for N = 20, spin coorindate axes 
for spheres as in Fig. [2] Nonadiabatic population loss P na 
from the cat state during the initial acceleration phase is 
shown magnified in the left inset (magenta). The back panel 
shows the interference signal found in the relative distance 
distribution p(z) at t = rtrap/2 (blue), note the different ab- 
scissa used. 



ro = 209.42/is [39]. For our parameters r\n e = 0.23ms, 
larger than the time required for initial state creation 
(0.15ms) and acceleration (6/is). 

Conclusion: We have proposed a setup in which two cold 
atom clouds of about 20 atoms each evolve dynamically 
by internal forces into a spatial Schrodinger's cat state 
if exposed to Rydberg dipole-dipole interactions through 
dressing. The interactions create a state where two en- 
tire atomic clouds simultaneously are at two quantum 
mechanically superimposed locations, which are macro- 
scopically distinguishable. Hence they can be resolved by 
visible light. 

The internal forces that induce motion of the atomic 
clouds are also instrumental in creating the required in- 
termediate hyperfine state | ^ rep ) . This state may have 
interesting applications by itself due to its entanglement 
structure between the two clouds. Finally, the hyper- 
fine state | \I/ r ep ) + | ^att ) prior to any spatial dynamics 
realizes a collective spin Schrodinger's cat state. 

We gladly acknowledge fruitful discussions with Igor 
Lesanovsky, Klaus M0lmer, Markus Miiller, Pierre Pillet, 
Thomas Pohl and Shannon Whitlock, and EU financial 
support received from the Marie Curie Initial Training 
Network (ITN) COHERENCE". 



SUPPLEMENTAL MATERIAL 

This supplemental material provides additional details 
regarding the definition and use of Born-Oppenheimer 
surfaces, the spin structure of initial states as well as our 
proposal for initial state creation. 

Motion on Born-Oppenheimer surfaces: We insert 
the expansion | \£ ) = 0k (R, i)\ k) into the time- 
dependent Schrodinger equation iHd\^f)/dt = H\*&) 
with the Hamiltonian H from Eq. (1) of the main article. 
Upon projection onto the electronic basis state | k ) this 
yields a system of coupled Schrodinger equations for the 
atomic motion in R and electronic dynamics 



ih- 



c (R,t) 



dt 



2M 



R 



V(R) 



t (R,t) 



£ £ A*(R)<k|*£*g|i>, (4) 
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The form Q is used in the main article. 

It is often also instructive to convert to a picture us- 
ing the Born-Oppenheimer separation. To this end we 
expand the total wave function as 



|*) = ^<A k (R,i)|^(R)>, 



(5) 



in terms of eigenstates of the electronic Hamiltonian 



# int (R)| w (R)) = [/,(R)lw(R)>- 



(6) 



Deriving the equation of motion from the Hamltonian us- 
ing this expansion leads to the Born-Oppenheimer sepa- 
rated version of the Schrodinger equation: 



2M 



V(R) + Z7,(R) 



i(R,t) 



^A Im (R)^ m (R,f), 



(7) 



where Aj m (R) are non-adiabatic coupling terms [30]. As 
long as these remain small, components of the many- 
body wave function with different / are effectively de- 
coupled. The eigenvalues Ui(H) form separate Born- 
Oppenheimer potential energy surfaces that are in our 
context most useful to anticipate the atomic dynamics: 

If the many-body wave function is localized initially 
near Rq in a narrow region |R — Rq| < cr of the 27V 
dimensional parameter-space, atom n will be accelerated 
along the downhill gradient i^ n (Ro) = — V rn £//(Ro) of 
the energy surface /. Calculating these gradients for 
the surfaces Z7 rep discussed in the main text, we find 
approximately F reP)n (R) « ±^N[dD(r) /dr] \ d , where 
the lower (upper) sign applies for an atom n in cloud A 
(B). Fat^n has reverse signs. These expressions become 
exact for small a and large N. The 7V 2 -scaling of U rep / att 
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reflects the number of interacting atom pairs, and the TV- 
scaling of F re p,n the number of atoms excerting a force 
on atom n. Importantly, the force in those two states 
will induce motion of the clouds as a whole, since it is 
almost equally strong for all atoms. 
Spin analogy: After elimination of the Rydberg states 
our atoms are described with just two essential states, 
hence a mapping to coupled spin- 1/2 particles is possible 
with ||) = | <7 ) , |t) = | ft ) . We can then define collective 
spin operators for each cloud S^ A ^ B ^ 

S (A/B) = ^ S ("), (8) 
neA/B 

where the individual spin operators s^ n ^ = 
(ai n \ ay U \ ai^ ) T act on atom n only, and o~j are 
Pauli matrices. 

For equal interactions between all pairs of atoms from 
different clouds, D n i = D, which is approximately re- 
alized due to a <C d, the interaction Hamiltonian i^ int 
takes the form 

Am = D(S { _ A) S {B) + S {A) S {B) ), (9) 

where S^ B ^ = si A ^ B ^ =L iSy A ^ B ^ are collective raising 
and lowering operators. Restricting ourselves to the same 
Hilbert space as in the main body of the paper (27V atoms 
in both clouds together, with equal numbers in | g ) and 
| ft)), we see that only collective spin states with total 
spin S = N/2 in A/B and magnetic quantum numbers 
mi^ = —mi B ^ = —N/2---N/2 have to be considered. 
In terms of the assignment of | ft ) states, = N^a — 

7V/2, where N^a is the number of atoms in state |ft) 
within cloud A. 

The analogy to a spin model allows a spin-sphere rep- 
resentation of hyperfine states in our system, for which 
we calculate the Q-mnction defined by Q(ip, 0)[^] = 
| ( ip, 6 | ^ ) | 2 for each surface element ((/?, 0) of the sphere. 
These are shown in Figs. 2 and 3. 

The underlying coherent spin states 



\<p,e) = (l + \ V \ 2 )-i J2 \ (s + m K +m \S,m), 

m=-S V ^ ' 

rj = - tan (0/2) exp (-iy) (10) 

are defined as common in coupled spin-systems (Ref. 

[26]). We used |S,m) = MT,N hA =m+N/2 I k >> where 
J\f normalizes the state. 

Chirped micro-wave pulse: In the main body of the 
text we have described how interaction of our system 
with a time-dependent micro- wave field, described by the 
Hamiltonian 

#ini(R-) = i^int(R-) + ^rf? 

H r{ = 5> rf (^V2 + h.c. - A rf (t)4£], (11) 

n 



can be used to adiabatically create the fully repulsive 
electronic state | ^f rep ) . As we are interested in fast 
micro-wave pulses to minimize the chance for sponta- 
neous decay and to avoid an onset of motion, it is crucial 
to know the energy gap between the state that is to be 
adiabatically followed and other states in the spectrum. 

The problem separates for the case of no interaction 
^int(R) = 0, and each atom can independently be in the 
two eigenstates of the field Hamiltonian, which in matrix 
representation reads 

We denote the eigenstates by | + ) or | — ) . These have 
energies E± = (— A r f ± + The gap between 

the highest energy state (all atoms in | + )) and the state 
adjacent in energy (one atom in | — ), rest in | + )) is thus 
AE = E+ - E_ = v/A r 2 f + Q 2 rf . 

For the opposite case without r.f. field, we can numer- 
ically solve i^i n t(R) for configurations as in Fig. 1 and 
find a gap g(N)D(d), where the prefactor g(N) is only 
weakly TV-dependent and approaches g(N) ~ 0.5 for N 
of the order 10 or larger. For small A/", e.g., g(4) = 1.2. 
Our analytical expression in Eq. 3 interpolates between 
these two limiting cases through simple addition, and is 
found to describe all inspected cases satisfactorily. 

While the final gap for Q r f — » is thus roughly inde- 
pendent of the number of atoms, the fidelity in numerical 
simulations nonetheless decreases slightly for larger A/", 
due to stronger non-adiabatic couplings. It could proba- 
bly be increased by pulse-shape optimization. 
Conditional phase flip: It was described in the main ar- 
ticle how the maximally repulsive hyperfine state | ^ rep ) 
can be obtained from a state of all atoms in | g ) using a 
chirped microwave field. This section supplies details on 
the subsequent step, the transfer to the state | ^ cat ) using 
the scheme of Ref. [32] . For this we assume that dressed 
dipole-dipole interactions are adiabatically removed, so 
that | \I/ re p ) is now given in terms of bare ground states 
| g ) and | ft ) . 

Let there be a control atom with two internal states 
|0)c, 1 1 )c? embedded in cloud A, which is otherwise 
unaffected by the creation of | ^ rep ) so that the process 
described above corresponds to | )c \ g ) — >■ | )c \ ^rep )• 
The control atom could be another species of atom, or 
one in two hyperfine states different from \g), \ ft ). 

Consider the following protocol: 

1. Apply a Rabi tt/2 pulse on the control atom to 
obtain (|0) c + |l)c)/>/2|tf rep / • 

2. Using the mesoscopic Rydberg quantum gate of 
Ref. [32], iff the control is in \l)c we now trans- 
fer all atoms in cloud A only from | ft ) into a third 
hyperfine state \h'). This creates (| )c\ ^rep ) + 

1 1 )c \ ^rep ))/V% where | ^rep ) is obtained from 
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I ^rep ) by replacing each atom in cloud A that was 
in state | h ) by one in state | h! ) . 

3. Apply a Rabi 2tt pulse between | h! ) and some 
auxiliary level to obtain a phase (—1) per atom 
that was in | h! ) . This has created the state 
(|0)c|*re P ) + |l)c|*^tt))/V2, as can be seen 
from Fig. 2 (a). 

4. Apply the gate again to return all atoms from | h! ) 
to | ft), resulting in ( | ) c I ^re P ) + 1 1 ) c I * att ) ) / V2 , 
with which we have reached our goal. 
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